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Summary 

A low Reynolds number k-e turbulence model and condi- 
tioned momentum, energy and turbulence equations were used 
to predict bypass transition heat transfer on a flat plate in a 
high-disturbance environment with zero pressure gradient. The 
use of conditioned equations has been demonstrated by other 
researchers to be an improvement over the use of the global- 
time-averaged equations for the calculation of velocity profiles 
and turbulence intensity profiles in the transition region of a 
boundary layer. The present work extends the approach of con- 
ditioned equations to include heat transfer and uses a modeling 
of transition events to predict transition onset and the extent 
of transition on a flat plate. These events, which describe the 
boundary layer at the leading edge, result in boundary-layer 
regions consisting of (1) the laminar, (2) pseudolaminar, 
(3) transitional, and (4) turbulent boundary layers. The modeled 
transition events were incorporated into the TEXSTAN two- 
dimensional boundary-layer code which is used to numerically 
predict the heat transfer. The numerical predictions in general 
compared well with the experimental data and revealed areas 
where additional experimental information is needed. 

Introduction 

The transition from laminar to turbulent boundary-layer flow 
effects great increases in the local wall shear stress and heat 
transfer. This effect of transition is especially critical for airfoil 
surfaces such as turbine blades where 50 percent or more of 
the blade surface can be in transition. The ability to predict 
the starting location and streamwise extent of transition is 
important to the determination of turbine blade heat transfer 
that critically affects longevity and engine performance. The 
accurate prediction of gas-side heat transfer on turbine blades 
will become more critical as the desired operating temperatures 
and stage loading levels of advanced turbine engines increase. 
The present methods for predicting transition shear stress and 
heat transfer on turbine blades are based on incomplete 
knowledge. Little is known about the transition process in an 
engine environment where disturbance levels are initially large. 

In such a large disturbance environment, traditional linear 
mechanisms are bypassed and finite, nonlinear effects must 
be considered. 

Modeling of this transition phenomena to predict wall shear 
stress and heat transfer on turbine blades must consider the 


effects of frec-stream turbulence, pressure gradient, streamwise 
curvature, surface roughness, wall and free-stream temperature 
ratio and flow disturbances (e.g., wakes). In the present work 
the effect of free-stream turbulence on the computation of 
transition heat transfer on a flat plate is investigated. 

Computational models for the calculation of wall shear stress 
and heat transfer in transition flows may be classified into four 
groups (Narasimha, 1985). These groups are linear 
combination models, algebraic models, differential models, 
and higher order models. Linear combination models calculate 
overall shear stress and heat transfer in terms of a linear 
combination of the intermittent laminar and turbulent 
properties of the transition boundary layer. Examples of this 
group are Dhawan and Narasimha (1958), Chen and Thyson 
(1971), and Dey and Narasimha (1988). Algebraic models use 
an algebraic model for the Reynolds stress of the time-averaged 
equations of motion. By gradually turning on the turbulent 
viscosity in proportion to the intermittency, it becomes possible 
to simulate the transition from a laminar to turbulent boundary- 
layer. Such a model requires knowledge of the start of 
transition, the transition length, and the transition path. 
Examples of this are Ccbcci and Smith (1974) and Gaugler 
(1985). Differential models describe the Reynolds stress of 
the time-averaged equations by the use of one- or tw'o-equation 
turbulence closure models. The one-equation model determines 
a turbulent velocity scale with a turbulent kinetic energy 
equation with the length scale determined algebraically. The 
two-equation model determines the length scale with an 
additional partial differential equation. The most well-known 
two-equation model is the Jones and Launder (1973) k-e model. 
Examples of this model for simulating transition in a high- 
disturbance environment are given by Rodi and Scheurer 
(1985) and Schmidt and Patankar (1988). Higher order models 
do not use the eddy viscosity concept to calculate Reynolds 
stresses (as required in the above three models), but rather 
use differential equations for each of the Reynolds stresses. 
Donaldson (1969) performed a two-dimensional boundary- 
layer analysis with higher order equations to determine the 
effect of large disturbances on transition. Differential models 
require a knowledge of the starting location of the initial 
profiles to determine the correct location of the beginning of 
transition. The work of Schmidt and Patankar (1988) is an 
example of the importance of locating these initial profiles in 
transition calculations. 

The present work uses a low' Reynolds number, k-e turbu- 
lence model for calculating the transition heat transfer in 
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environments of high disturbance. In general, k-e two equation 
models simulate the transition governed by the transport of 
turbulence from the free stream into the boundary layer. 
However, previous results of this approach foi flat-plate heat 
transfer have underpredicted the transition length and incon- 
sistently predicted the start and finish of transition (Rodi and 
Scheurer, 1985). To address these problems and to make pre- 
dictions consistent with the experimental evidence, Schmidt 
and Patankar (1988) modified the turbulent production term 
in the differential equation for the turbulent kinetic energy (k) 
and assumed the presence of an initial laminar boundary layer. 

The present approach divides transition into four regions 
for the purpose of applying a model to the physics of each 
region: (1) the laminar boundary layer (2) the pseudolaminar 
(3) the transition and (4) the turbulent boundary layers. The 
conditions for determining the existence of each region will 
be presented. Models will be given for determining disturbance 
growth in region 2 and describing the presence of turbulent 
spots in region 3. Numerical predictions are made using the 
TEXSTAN boundary-layer computer code (Crawford, 1985). 
TEXSTAN is based on the STAN5 boundary-layer program 
developed by Crawford and Kays (1976). The finite difference 
scheme of TEXSTAN is based on the numerical algorithm by 
Patankar and Spalding (1970). TEXSTAN solves the steady 
two-dimensional parabolic differential equations that govern 
boundary-layer flow. This program sequentially solves the 
momentum equation and any number of transport equations such 
as stagnation enthalpy, turbulent kinetic energy (TKE), turbulent 
dissipation rate (TDR), and mass concentration governing 
equations. Numerical predictions are compared with the flat* 
plate, zero pressure gradient data of Blair and Werle (1980). 

Symbols 

C /<( C|,C 2 constants appearing in the k-e turbulence model 
model constants for source terms 
skin friction 

low Reynolds number functions 
intermittency function 
thermal conductivity 
turbulence kinetic energy 
dissipation length scale 
transition length 

nondimensional spot formation rate 
spot formation rate 
turbulent Prandtl number 
static pressure 
wall heat flux 

momentum thickness Reynolds number 
Stanton number 
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/ 

K 

k 

4 

Ur 
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n 

Pr, 

P 

Rc 0 

St 


temperature 
fluctuating temperature 
turbulence intensity 
turbulent heat flux 
time 

free*stream velocity 
mean velocity in x,y,z directions 
' t w' fluctuating velocities in x,y,z directions 
turbulent shear stress 

location of transition onset (spot formation begins) 
y intermittency 

6* displacement thickness 

§ 3 energy thickness 

e dissipation rate 

S boundary-layer momentum thickness 

A longitudinal integral length scale 

\ thermal conductivity/specific heat 

ix molecular viscosity 


T 

T 

Tu 

7V 

l 

U c 

U,V,W, 
ll' ,Y 
14 ' V ' 
-t/r 


Pt 

eddy or turbulent viscosity 

p 

fluid density 

a 

dependence area factor 


empirical constants in turbulence n 

Subscripts: 


e 

denotes free-stream value 

l 

laminar zone of transition region 

t 

turbulent zone of transition region 

Jr 

transition 


Analysis 

The boundary layer flow is assumed to be steady, incom- 
pressible, and two-dimensional. Spanwise variations are 
neglected. 

Figure 1 depicts the development of a boundary layer in x-y 
coordinates. For purposes of modeling, the growth of the 
boundary layer is described in terms of regions. Region 1 is 
the laminar boundary layer, where small disturbances intro- 
duced into this layer are damped. Experimental evidence 
suggests that this region is nonexistent for free-stream turbu- 
lence levels greater than 5 percent. Region 2, which begins 
when the critical Reynolds number for linear instability is 
reached, represents the growth of nonlinear disturbances and 
extends from the critical Reynolds number to the position where 
the turbulent spots begin. In region 2, according to Morkovin 
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Region 

1 Laminar boundary layer 

2 Psuedo-Laminar 

3 Transition 

4 Turbulent boundary layer 



y=o 

Figure I —Transition model. 

(1978), when large disturbances exist, linear stability mecha- 
nisms are bypassed and finite, nonlinear instabilities occur. 
Because this region may have characteristics of a laminar bound- 
ary layer, it is called the pseudolaminar region. The growth 
of unstable waves in region 2 eventually results in the formation 
of turbulent spots (region 3). In region 3, the beginning and devel- 
opment of turbulent spots, there is an intermittent appearance 
of turbulent spots that grow as they move downstream until they 
finally merge to form the turbulent boundary layer (region 4). 
Each of these regions will be considered separately. 

Region 1 (Laminar Boundary Layer) 

The work of Abu-Ghannam and Shaw (1980) indicates that 
the transition momentum Reynolds number decreases to an 
asymptotic value of 163 for increases in the free-stream 
turbulence intensity. This Reynolds number value suggests that 
region 1 is governed by considerations of linear stability or 
Tollmien-Schlichting stability. Contrary evidence for high 
disturbance levels is indicated by some recent work (Rued and 
Wittig, 1985; Blair and Werle, 1980; Sohn, Reshotko, and 
O’Brien, 1989). These results indicate that the onset of 
transition begins close to the leading edge of a flat plate for 
free-stream turbulence intensity greater than about 5 percent. 
Rued and Wittig theorized that the lack of a minimum 
asymptotic transition Reynolds number for his results was due 
to the use of a thin leading edge. They speculated that the 
acceleration effects of thick leading edges provide additional 
stability. Such differences would be important when applying 
flat-plate transition information to turbine design. 

Elder (1960) demonstrated that below the critical Reynolds 
number for stability, the flow was stable to small disturbances 
but unstable to large disturbances. He determined that 
breakdown to turbulent spots would occur for boundary-layer 
velocity fluctuations greater than 18 percent of the free-stream 
velocity. He found this event to be independent of the Reynolds 
number, and he postulated that the Reynolds number was 
important in the growth of the disturbances. 


Dyban’s experiments (1976) on a flat plate suggest that 
disturbances of the order determined by Elder for transition 
to occur are present near the leading edge for free-stream 
turbulence intensities of 4V2 percent or more. The evidence, 
therefore, suggests that in cases where additional stabilization 
(e.g., acceleration) is not present, region l is essentially 
nonexistent for free-stream intensities of 5 percent or more. 
This is consistent with the statement of Reshotko (1986) that 
transition can occur at Reynolds numbers lower than the 
Reynolds for the onset of linear growth when initial disturbance 
is large enough to be nonlinear. According to Reshotko (1986) 
who cites the experimental data of Sohn (1986), the onset of 
nonlinear boundary-layer disturbance growth for 2 percent 
free-stream turbulence begins at a momentum Reynolds 
number of 96. This suggests that the onset of nonlinear grow th 
is dependent on the level of the free-stream turbulence. 
However, the transition Reynolds number of 96 measured by 
Sohn is too low. Sohn attributes his value to possible adverse 
pressure gradient effects. Until definitive experimental 
evidence is provided, the present work assumes that nonlinear 
boundary-layer disturbance growth begins at the Reynolds 
number for the linear stability limit. 

Critical momentum Reynolds numbers of 163 and 200 are 
generally used for the onset of instability, with 200 being 
considered the more accurate value. However, for the present 
study a value of 163 is used for consistency with the asymptotic 
limit at high turbulence intensities reported by Abu-Ghannam 
and Shaw (1980). 

The following criteria are used for region 1 : 

1. No stabilizing effects present at the leading edge 

a. Tu (i < 5 percent; disturbance growth (region 2) 
begins at Re 0 - 163 

b. Tu e > 5 percent; disturbance growth begins at 
leading edge 

2. Stabilizing effects present at leading edge; disturbance 
growth begins at Re tf =163. 

Region 2 (Pseudolaminar) 

Region 2 requires that we analytically describe the 
amplification of boundary-layer disturbances in conjunction 
with the momentum and energy equations for the boundary 
layer. Dey and Narasimha (1984) proposed that transition is 
turbulence-driven for turbulence intensity levels greater than 
0. 1 percent. For levels less than this value, transition is said 
to be controlled by the availability of nonturbulent disturbances 
such as noise and vibration. This suggests that bypass transition 
and the amplification of disturbance energies (Reynolds 
stresses) can be described by turbulence models which 
determine the convection, diffusion and production of 
turbulence energy in the boundary layer. The equations to be 
used for region 2 are the boundary-layer equations for 
momentum and energy and the Jones-Launder (1973) two- 
equation turbulence model. These equations are given in the 
following sections: 
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Boundary Aayer equation .— The flow is assumed to be 
steady and two-dimensional neglecting spanwise variations. 


TABLE I.— CONSTANTS AND FUNCTIONS FOR JONES-LAUNDER 
TURBULENCE MODEL 


Momentum equation: 


du du d ( du —— \ dp 

pu 7" + pv V = T t V 7" “ pU v ) ~T 

ox ay dy \ oy / dx 


( 1 ) 


where 




by \ e J by 


( 2 ) 


Energy equation: 


dT dT d t dT -y—, 

pu 1 - pv — = — A pv T 

dx dy dy V dy 


(3) 


where 


-rzr, Pi dT 

-pv T = 

Pr, dy 


(4) 


Jones-Launder turbulence model — Turbulence kinetic 
energy: 


dk dk d 

pu I- pv — -= — 

dx dy dy 


" oj dy 


>Y 

■ * w 


- pe - Ip 


dk 


l/ 2\2 


dy 


(5) 


Turbulence dissipation: 


de de d 

pu 1- pv — — — 

dx dy dy 


, Pt\ 

M oj dy 


e f du\ 

+ V'U. 


Cif^pz 2 - „ ppi ( drii\ 

— - - + 2.0 — ( — t I (6) 

k p \dy 2 J 


Constants for the turbulence model are give in table I. 

Region 3 (Transition) 

Region 3, the transition region, is characterized by the 
intermittent appearance of turbulent spots which grow as they 
move downstream until they finally merge to form the 
turbulent boundary layer. The beginning of the transition 
region is defined by the first appearance of turbulence spots. 
Numerically this is assumed to occur when the calculations 
for region 2 show' a rapid increase in turbulence kinetic energy 
k which results in an increase in the skin friction. Vancoillie 
and Dick (1988) state that conventional turbulence models 


(a) k-c Turbulence model constants 


C M 0.09 

C, 1.45 

C 2 2.0 

a, 10 


(b) Functions 

/« 

exp [-2. 5/(1 + /?,/50)] 

f 2 

1 .0 - 0.3 exp - R 2 j 

* T 

pitl/p. = pk 2 /ep 


based on global time averages can not give a good description 
of such intermittent flow, resulting in poor agreements for the 
turbulence intensity profiles for the boundary layer. They 
reported that using conditional averaging techniques for the 
turbulent spots and the laminar-like fluid surrounding the spots 
resulted in a good prediction of the experimental values of the 
velocity profiles and the turbulence intensity profiles in the 
transition region. As stated previously, turbulence models 
based on global averages underpredict the transition length. 
Schmidt and Patankar (1988) corrected this by using an 
empirical modification of the turbulence production term. 

The present work extends the conditioned equation method 
of Vancoillie and Dick to include the energy equation for 
predicting heat transfer in the transition region. The condi- 
tioned momentum, energy, and turbulence equations for the 
turbulent spot and nonturbulent (laminar) portion of the 
intermittent flow in the transition region are combined and 
simplified to obtain global values of velocity and temperature. 
This permits use of the boundary-layer TEXSTAN code. 
Details of the analysis are presented Tn the appendix. The 
equations developed in the appendix (eqs. (A- 12) and (A-23)) 
permit a consideration of the Reynolds stresses in the 
nonturbulent (laminar) zone. However, as assumed by 
Vancoillie and Dick and verified by the works of Sohn and 
Reshotko (1991) and Kim, Simon, and Kestoras (1989) for 
a flat plate with zero pressure gradient, the Reynolds stresses 
of the laminar zone can be neglected. As shown in the 
appendix, when the modeled equations for the turbulent and 
laminar zones are combined and the assumption made of 
negligible Reynolds stresses in the laminar zone, the following 
momentum and energy equations are obtained: 


du du 

pu F pv — - 

dx dy 


d 

d: V 



- 7P(«A V) 


dp 

dx 


+ 7 p(w/ - «/) 
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dx 


du 

dx 
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dy 
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dT 
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d dT 
dy dy 
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Pr, dy 
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dT 
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dT 
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(8) 


The Jones-Launder turbulence model, as derived in the 
appendix (eqs. (A-26) and (A-27)) for the zone of turbulent 
spots, in terms of global averages for velocity, is as follows 
with the assumptions of equation (9): 


Equation (7) shows the Reynolds stresses multiplied by the 
intermittency factor. This is a form used by McDonald and 
Fish (1973) and others. A comparison of equations (7) and 
(I) indicates additional source terms as a result of spot 
formation. A comparison of equations (8) and (3) indicates 
the need to consider additional source terms because of the 
presence of a turbulent zone and a laminar zone in the transition 
region. These equations are simplified to permit the use the 
boundary-layer code TEXSTAN. The simplified equations also 
permit an evaluation of the additional source terms and of the 
validity of simplifying assumptions. The following assumptions 
are made for the momentum and energy equations: 

du, _ du du, du 

dy dy dx dx 

(9) 

dT } _ dT dT, _ dT 

dy dy dx dx 

The assumption for the velocity profile was also made by 
McDonald and Fish (1973) in their development of a turbu- 
lence model for transition. This assumption becomes more 
applicable as intermittency increases and velocity profiles 
become more turbulent. The errors introduced by this assump- 
tion should not be too large in the early part of transition 
because of the small value of intermittency. The result of these 
assumptions, as applied to equations (7) and (8), is the 
modification of equations (1) and (3) with the intermittency 
factor as follows: 
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(13) 


Equations (10) to (13) are the simplified equations for the 
transition region. These equations require the specification of 
intermittency. 

Specification of intermittency requires knowledge of the 
transition path in terms of the transition start and length. 
Narasimha (1957) derived from the turbulent spot theory of 
Emmons (1951) the following transition path equation: 


7 = 1.0- exp 



for 7 = 0 — 0.99 


(14) 


Using the approach of Narasimha (1985) the transition length 
may be expressed in terms of the transition Reynolds number 
and a nondimensional spot formation rate as follows: 


2.15 

Re, =— Re« /2 


v9v 


-v lr 


(15) 


where N, the nondimensional spot formation rate is defined as 


where 


N = nodlrp/p 


(16) 


du 


-p(w/v/), = #*,— = 

dy 


( C^k} \ du 

\ * )dy 


and 


Narasimha (1985) found that for free-stream turbulence 
levels greater than 0. 1 percent the value of vV has the approxi- 
mate constant value of 0.7 x 10 ~ 3 . Equation (15) was used 
for determining the value of N for free-stream turbulence levels 
greater than 0.3 percent. For the data of Sohn and Reshotko 



0 2 4 6 8 10 12 14 

x, in. 


(a) Free-stream turbulence intensity, Tu,„ 1.1 percent; transition length, 
L tr 1.03 ft; nondimensional spot formation rate, 0.53 x 10 momentum 
thickness Reynolds number, Re 0 , n 337. 

(b) Free-stream turbulence intensity, Tu r , 2.4 percent; transition length, L ir 
0.81 ft, nondimensional spot formation rate, 0.26 x 10 \ momentum 
thickness Reynolds number. Re#,,., 229. 

Figure 2. — Narashimha method for determining 7 = 0 and 0.99 points. (Data 
from Sohn and Rcshotko (1991).) 



Free-stream turbulence intensity, Tu e , percent 

Figure 3. — Variation of nondimensional spot formation rate with free-stream 
turbulence (Narasimha (1985) method). 


(1991), the transition start and length were determined by 
plotting the intermittency in the manner of Narasimha (1957) 
and determining the distances x for the 0 and 0.99 intermittency 
points (fig. 2). The data of Sohn and Reshotko are for cross- 
stream measurements of turbulence. It is assumed that the 
turbulence of their experiments was nearly isotropic such that 
their measured turbulence differs little from the total turbulence 
as defined in this report. The present work uses a lower value 
for N of 0.3 x 10 ” 3 , as indicated in Figure 3. The use of this 
value in equation (16) results in the following equation for the 
calculation of the transition length: 

Re L = 124 Re* 3/2 for y = 0 — 0.99 (17) 

hr u rr 

Region 4 (Turbulent Boundary Layer) 

Calculations for this region are made using equations (11) 
and (12) with the intermittency set at one. 

Initial and boundary conditions. — Wall boundary condition: 
Along the wall the no-slip boundary condition was applied as 
follows: 

j = 0:m = v = 0 (18) 

For the energy equation a wall heat flux was prescribed. 

q w =-lP(y = 0) (19) 

dy 

The turbulence kinetic energy and dissipation rate are set 
to zero at the wall. 

y = 0 : * = e = 0 (20) 

The zero boundary condition for the dissipation was made 
possible by Jones and Launder (1970) who added terms to the 
turbulence energy equation (eqs. (5) and (12)). 

Free-stream boundary conditions: The velocity (£/,,) and 
temperature ( T t ) were set equal to a constant along the free 
stream. Boundary conditions for the turbulence kinetic energy 
and dissipation rate are determined from solving the transport 
equations (5) and (6) at the edge of the boundary layer. 


dk c _ 

e<< 

(21) 




(22) 

dx k„ 



where initial values of the free-stream turbulence kinetic 
energy are determined from the turbulence intensity as follows: 

k e = 1.5(Tu f I/ f ) 2 (23) 
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where 



and the initial dissipation rate is determined from 



c 

The dissipation length scale L e can be calculated from the 
longitudinal integral length scale. 

Initial values.— To initiate calculations in region 1, a Blasius 
velocity profile and a flat temperature profile (unheated start) 
were used. For the turbulence transport equations, the 
following initial profiles suggested by Rodi and Scheuerer 
(1985) were used: 

k = k c (u/U c ) 2 (25) 

and 


e = 0.35k — (26) 

dy 

Results and Discussion 

Calculations were performed using the TEXSTAN code for 
turbulence levels of 1 .4, 2.8, and 6.2 percent. The calculations 
were chosen for comparison with the experimental data of Blair 
and Werle (1980). They conducted their experiments in a low- 
speed wind tunnel (100 ft/sec) under ambient air conditions 
and used grids to generate turbulence. The test section was 
heated at a rate of 0.078 Btu/ft 2 - sec w'ith an unheated length 
of 0. 141 ft. The experimental boundary conditions used in the 
calculations are listed in table II. According to Hancock 
(1980), for grid generated turbulence, the measured values 
of longitudinal integral length scale can be used to calculate 
the dissipation length scale as follows: 

L c = 1.5 A, (27) 


TABLE tl. -BOUNDARY CONDITIONS FOR BLAIR- WFRLE DATA 


Parameter 

Grid I 

Grid 2 

Grid 3 

Free-stream velocity, U r , ft/s 

100.0 

100.0 

100.0 

Free-stream temperature, T ( ., "F 

68.5 

68.5 

68.5 

Wall heat flux, r/ w , Blu/fr/s 

0.078 

0.078 

0.078 

Free-stream turbulence intensity, Tu^, percent 

1.4 

2.8 

6.2 

Free-stream turbulence kinetic energy, k <Mt ft"/ s 2 

2.94 

11.8 

57.7 

Free-stream dissupation rate, c ( . y ft Vs 1 

127 

511 

2450 



0 1 2 3 4 5 6 7 

Free-stream turbulence intensity, Tu e , percent 

Figure 4.— Comparison of predicted transition onset with experiments. 


Transition Onset 

The results of the calculation for the beginning of transition 
using the equations for region 2 and the criteria for the start 
of region 2 are given in figure 4. As indicated previously, the 
initiation of turbulent spots (beginning of region 3) is assumed 
to occur when the calculations for region 3 show a rapid 
increase in the skin friction. Also shown in figure 4 are the 
experimental results of Gardiner (1987), Gostelow and Blunder 
(1988), Sohn, Reshotko, and O'Brien (1989) and the experi- 
mental correlation of Abu-Ghannam and Shaw (1980). 
Figure 4 shows the effect of assuming stabilizing effects at 
the leading edge for turbulence levels greater than 5 percent. 
For the 6.2-percent turbulence level, the present calculations 
agree with the correlation of Abu-Ghannam and Shaw, if we 
assume stabilizing influences at the flat-plate leading edge. If 
we assume that boundary-layer instability begins at the leading 
edge, and not at the critical Reynolds number for instability, 
the result is the lower value for transition indicated in figure 4. 
This is in good agreement with the measured results of 
Gostelow and Blunder and the observations of Rued and Wittig 
(1985), Blair and Werle (1980), and Sohn, Reshotko, and 
O’Brien (1989). As indicated previously, the criterion for the 
instability beginning at the leading edge for free-stream 
turbulence levels greater than 5 percent was used in the present 
calculations. Figure 4 shows general good agreement of 
experiment and calculations. Gardiner (1987) and Gostelow 
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Stanton number, St 


.0040 | — 



(a) Total turbulence intensity, Tu t ,, 1.4 percent. 

(b) Total turbulence intensity, Tu <( , 2.8 percent. 

(c) Total turbulence intensity, Tu ( „ 6.2 percent. 

Figure 5. —Comparison of prediction with experiment (zero pressure gradient). 



Figure 6.— Skin friction distribution for total turbulence intensity of 6.2 percent. 


and Blunder (1988) defined the transition at the 1 percent 
intermittency level while the transition values of Sohn and 
Reshotko (1991) are for zero percent intermittency. According 
to the numerical calculations, this difference in definition 
should result in a small difference in the Reynolds number 
(approximately 10 percent), as confirmed by the slightly lower 
experimental values of Sohn and Reshotko as compared with 
those of Gostelow and Blunder. As expected, the experimental 
correlation of Abu-Ghannam and Shaw indicates higher transi- 
tion values than other experiments, since their definition for 
transition was based on the onset of deviation from laminar 
boundary-layer characteristics. 

Transition Region 

Using the present analysis, as described above, for the begin- 
ning of transition permits a determination of when to begin 
the calculation for region 3. The resulting calculations, using 
equations (7) to (17), are compared with Blair and Werle’s 
experimental data in figure 5. There is generally good 
agreement, although the calculated curve swings below the 
experimental data. This may be the result of neglecting the 
source terms in the momentum and energy equations (eqs. (7) 
and (8)) and a result of the presence of high levels of low- 
frequency unsteadiness in the nonturbulent (laminar) zone of 
the intermittent boundary layer (Sohn, Reshotko, and O’Brien, 
1989), which the models of the present investigation may not 
have taken into account. Even after spot formation begins, the 
boundary layer acts as if it were a laminar boundary layer up 
to intermittency values of 0.40 (1.4 percent case) and 0.49 
(2.8 percent case). This finding is consistent with the measured 
velocity profiles of Sohn, Reshotko, and O Brien (1989), 
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which showed a Blasius profile for the laminar zone and a 
laminar-like profile for the overall profile, for intermittency 
values up to 0.34 with a 1 percent free-stream turbulence. This 
finding is also consistent with Gaugler (1985), who had to force 
transition early in what appeared to be the laminar region in 
order to agree with the experimental data. With intermittency 
values greater than 0.34, there is increased deviation from the 
Blasius velocity profile in the laminar zone. For this reason 
Sohn, Reshotko, and O Brien (1989) refer to the laminar zone 
as a nonturbulent zone to distinguish it from a purely laminar 
zone. 

The increase in heat transfer at substantial values of the 
intermittency highlights the importance of defining the 
transition point, since Blair and Werle defined transition as 
the location where the Stanton number first exceeded the 
laminar heat-transfer rate. The calculated transition point for 
the turbulence level of 6.2 percent (fig. 5(b)) occurs before 
the start of heating; therefore, there is no minimum in the 
Stanton number. This minimum becomes apparent if the 
calculated skin friction in the unheated length is plotted (fig. 6). 


Concluding Remarks 

The present work has examined the ability of two-equation 
turbulence models to simulate the onset of bypass transition 
and the region of transition consisting of islands of turbulent 
spots in a sea of laminar-like fluid. A critical test for a 
turbulence model is its ability to simulate disturbance growth, 
which is governed by the free-stream turbulence energy, and 
to determine the onset of transition. The Jones-Launder 
turbulence model has performed this job well and suggests that 
models for bypass transition need to comprehend the 
relationship between disturbance growth and the convection, 
diffusion, and production of turbulence in the boundary layer. 
The intermittent flow of the transition region was modeled by 
the use of conditioned flow equations rather than by global 
time-averaged equations. These equations better express the 
physics of turbulent spot-fluid interaction and produce a more 
complete and accurate description of turbulent energy budgets. 
In addition, they bring out the necessary and strong role that 
intermittency plays in calculating heat transfer with transition 
without the need to modify the turbulence production term as 
others have done (Schmidt and Patankar (1988)). 

It appears that the existence and length of the laminar 
boundary layer at the leading edge of a Hat plate depends on 
the presence of stabilizing or destabilizing effects. This 
suggests that application of the present work to turbine blades 
will require additional information on and modeling of the 
effect of curvature and pressure gradient on the stability of 
the laminar boundary layer. In the present effort for flow over 
a flat plate, it is assumed that nonlinear growth for free-stream 
turbulence levels less than 5 percent begins at the critical 
Reynolds number for linear stability. It is not clear that one 
is justified in this assumption. It is quite possible, as previously 


suggested, that the condition for the beginning of nonlinear 
growth is dependent on the free-stream turbulence level in 
addition to the effects of curvature and pressure gradient. There 
is a need in this area for more experimental and analytical 
information. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 28, 1991 

Appendix— Transition Region Equations 

The development of equations for the transition region uses 
the approach of Vancoillie and Dick (1988), who developed 
conditioned continuity, momentum, and turbulence equations 
for the intermittent flow in the transition region. In this 
appendix the conditioned equations for the energy equation 
are derived and the conditioned equations are combined to 
permit the calculation of global values. 

Vancoillie and Dick used an intermittency function I(x,y,z,t) 
for their conditioning analysis, with the function having a value 
of 1 inside the turbulent spot and 0 outside. The time average 
of I is the intermittency factor. Conditioned averages are 
defined by 


— ul — w(l -/) 

w, = - W/ = — (A- 1) 

7 1-7 

where 


7 = “ Idt = / (A-2) 

/ Jo 

By using equation (A-l), the following global averages may 
be developed: 


u = yu t + (1 - y )Uf 


(A-3) 


and 


v “ 7 V / + (1 -7 ) v, (A-4) 

In what follows, the use of a bar for simple symbols is omitted 
for convenience. 


Continuity Equations 

Vancoillie and Dick’s derivation of the continuity equations is 


dyu { + dyv t dl 
dx dy dt 


(A— 5) 
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result in the equation of global -averaged and conditioned- 
averaged variables: 


d( I - y)u, + <>(\ - y)vi 
dx ■ dy 


ill 

It 


(A-6) 


du dit 

' , "S + ' n '5T 


dp d 
— + — 
dx dy 



7 P(wA’/)/ 


Equations (A-5) and (A-6) may be combined, and using 
equation (A-2) and (A-3) results in the following global 
average continuity equation: 


- (1 - y)p(ufvf )/ 


+ pi 


(hi - ‘I,) 




t + t = 0 

dx dy 


(A-7) 
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( du, 

du\ 

W 
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(A- 12) 


Momentum Equations 

The modeled momentum equations of Vancoillic and Dick are 


du, 

py*<i — + py y * 
dx 


du, 

dy 


dp d 

-y — + 7 — 
dx dx 


p p(w f v, )/ 

av 


dy 

+ p(W/ - W f ) — + Cyfl 

dx 


u { - u, dy 
6* dx 


(A-8) 


Energy equations 

The procedure for determining the time averaged condi- 
tioned energy equations follows that used by Vancoillie and 
Dick for the conditioned momentum equations. The Navier- 
Stokes energy equation is multiplied by I (turbulent zone) and 
l - / (laminar zone), and the total derivative for I is multiplied 
by 77 and 7/ for the turbulent and laminar zone, respectively. 
The equations are added and time averaged. Some terms are 
neglected on the basis of boundary-layer approximations and 
two-dimensional flow. The result is as follows: 


and 


dui 


duj 

p( 1 - 7 )«/ 7T + “ 7) 1 
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dx 
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ay 
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Uf — u, dy 
b* dx 


dT t dT, 

P7»f — + p y Vf !T 
dx dy 




P(v/7T), 


3773 / 37/3/ T ,<// 

X — - — - A + pT t — 

dx dx dz dz dt 


(A- 13) 


for the turbulent zone and 


(A-9) 


Vancoillic and Dick solved equations (A-8) and (A 9) 
(neglecting the Reynolds stresses of laminar zone) to predict 
for the transition region, boundary-layer thickness factors, 
velocity, and turbulence intensity profiles. Since it is the objec- 
tive of the present work to determine the combined effect of 
turbulent spots and the laminar zone on heat transfer in the 
transition region, the individual equations for each zone are 
combined to obtain global values. 

Combining equations (A-8) and (A-9) and using equations 
(A-3) and (A-4) and the following derivatives of equations 
(A-3) and (A-4) 


dTt 37) /f 

p(l - y)ui — + p(l - y) v t — 0 

dr oy 



X 


37) 

3v 


1 37 / 3 / , 37/3/ di 

T!) - +x j^ +x T : »r eTi 7 , 


(A- 14) 

for the laminar zone. 

The first two source terms in equations (A- 1 3) and (A- 14) 
may be modeled 'in the manner of Vancoillie and Dick, 
resulting in 


du du, duf dy 

— - y — + (1 - 7) — + ( w / - M /) ~T 

dx dx dx dx 

du du, du, 

— = 7 — + (> _ T) — 
dy dy dy 


(A- 10) 
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and 


result in an energy equation expressed in terms of global and 
conditional averages as follows: 


d T, 

p( I - 7 )«/ — + p( 1 - y)v/ 


d.x 


dT, 

ay 


= d 



A^-p(v-/7/), 

dy 



T.) 


dy 

dx 


+ pli- 

ch 


(A- 16) 


Combining equations (A- 1 5) and (A- 16) results in 
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Jones-Launder Turbulence Model 

The Jones-Launder turbulence kinetic energy model derived 
by Vancoillie and Dick for the turbulent zone is 



P(v,'77), 


+ 




P ( v( 7/ ), | 


„ <// dl 

pt;-- p t;- (a- 17) 

di clt 


The source terms at the end of equation (A- 1 7) result from 
the interaction between the turbulent and laminar zones. These 
terms are modeled according to the method of Vancoillie and 
Dick. They assumed that the intermittency function at the 
interface between the laminar zone and turbulent spot changed 
linearly from 0 to 1 over a finite interval that goes mathe- 
matically to zero and that the flow velocity at the interface 
is equal to the laminar zone averaged velocity Uj. With these 
assumptions the above source terms may be expressed as 
follows: 


dl dl dy 

pT<- ~ pTf — = (7, - 7>, -f 
dt dt dx 


(A- 18) 


Treating equation (A- 17) in the same manner as the momentum 
equations (A-8) and (A-9) and using the following equations 

T — yT, + (1 — y)T, (A- 19) 
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It follows that the equation for the laminar zone is 
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Equations (A-24) and (A-25) are combined and simplified as 
previously. To permit simplification, we assume 

dk, _ dk, 

Av ~ 3y 

dk,_dk, 

dy dy 


li 


This suggests that the disturbance energies in the laminar zone 
are related linearly to those occurring in the turbulent spots. 
Experimental data (Sohn, Reshatko, and O'Brien, 1989, and 
others) seem to confirm this. The turbulence dissipation and 
turbulence viscosity terms for the laminar zone are small 
compared with the same terms for the turbulent zone and are 
therefore neglected. The result of combining is 
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Approaching the turbulence dissipation equations in the same 
manner as for the turbulence kinetic energy we have 
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